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Abstract 

A flow of metrics, gt, on a manifold is a solution of a differential equation d t g = S(g) , 
where a geometric functional S(g) is a symmetric (0, 2)-tensor usually related to some kind 
of curvature. The mixed sectional curvature of a foliated manifold regulates the deviation of 
leaves along the leaf geodesies. We introduce and study the flow of metrics on a foliation, 
whose velocity along the orthogonal distribution is proportional to the mixed scalar curvature, 
Sc m i x . The flow is used to examine the question: Which foliations admit a metric with a given 
property of Sc m j X (e.g., positive or negative)? We observe that the flow preserves harmonicity of 
foliations and yields the Burgers type equation along the leaves for the mean curvature vector H 
of orthogonal distribution. If H is leaf-wise conservative, then its potential obeys the non-linear 
Schrodinger heat equation. Under certain conditions the asymptotic behavior of its solutions is 
the same as of the linear equation, when the Burgers equation for H has a single-point global 
attractor. The flow admits a unique global solution g t , whose Sc m i x converges exponentially to 
a leaf- wise constant. The metrics are smooth on M when all leaves arc compact and have finite 
holonomy group. Using rescaling along D one may provide convergence to metrics with Sc m ; x 
as positive so negative. 

Keywords: foliation; flow of metrics; conformal; mixed scalar curvature; mean curvature; 
Burgers equation; non-linear Schrodinger heat equation. 
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Introduction 

1. Geometry of foliations. Let M n+P be a connected closed manifold (i.e., compact without 
a boundary), endowed with a p-dimensional foliation J 7 , i.e., a partition into submanifolds (called 
leaves) of the same dimension p. We have the g-orthogonal decomposition TM = Djr D, where 
the distribution Djr is tangent to T . Denote by ( • ) JF an d ( ■ )" L projections onto Dj- and D, 
respectively. The second fundamental tensor and the mean curvature vector field of J- are given by 

h T (X,Y) :=(l/2)(V x L' + VyV) ± , H T = Tr g h (X,Y e D T ), 

where V is the Levi-Civita connection of g. A Riemannian manifold (M,g) may admit many 
kinds of geometrically interesting foliations. Harmonic and totally geodesic foliations are among 
these kinds that enjoyed a lot of interest and investigation of many geometers (see survey in [8]). 
A foliation with vanishing mean curvature of the leaves is called harmonic. A foliation is totally 
geodesic if the leaves are totally geodesic submanifolds (i.e., the second fundamental form of the 
leaves vanishes) . The simple examples of geodesic foliations are parallel circles or winding lines on 
a flat torus, and a Hopf field of great circles on the sphere S 3 . The second fundamental tensor and 
the integrability tensor of D are given by 

h(X,Y) := (1/2) (VxY + VyXf, T(X, Y) := (1/2) [X, Yf (X,Y € D). (1) 
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The mean curvature vector of D is given by H = Tr g h. A foliation T is said to be transversely 
harmonic, or Riemannian, if, respectively, H = 0, or h = 0. Conformal foliations, i.e., h = 
(l/n)H ■ <7m, were introduced by Vaisman [T3j as foliations admitting a transversal conformal 
structure. Such foliations extend the class of Riemannian foliations. 

One of the principal problems of geometry of foliations reads as follows, see 
Given a foliation T on a manifold M and a geometric property (P), does there exist a Rieman- 
nian metric g on M such that T enjoys (P) with respect to g? 

Such problems of the existence and classification of metrics on foliations (first posed explicitly by 
H. Gluck for geodesic foliations) have been studied intensively by many geometers in the 1970's. 
A foliation is geometrically taut, if there is a Riemannian metric making T harmonic. H. Rummler 
characterized such foliations by existence of an ^-closed p-form on M that is transverse to T . 
D. Sullivan provided a topological tautness condition for for geometric tautness. By the Novikov 
Theorem (see [3]) and Sullivan's results, the sphere has no 2-dimensional taut foliations. In recent 
decades, several tools for proving results of this sort have been developed. Among them, one may 
find Sullivan's foliated cycles and new integral formulae, see [11[ Part 1], [15] etc. 

2. The mixed scalar curvature. There are three kinds of Riemannian curvature for a folia- 
tion: tangential, transversal, and mixed (a plane that contains a tangent vector to the foliation and 
a vector orthogonal to it is said to be mixed). The geometrical sense of the mixed curvature follows 
from the fact that for a totally geodesic foliation, certain components of the curvature tensor, see 
[8j. regulate the deviation of leaves along the leaf geodesies. In general relativity, the geodesic devi- 
ation equation is an equation involving the Riemann curvature tensor, which measures the change 
in separation of neighboring geodesies or, equivalently, the tidal force experienced by a rigid body 
moving along a geodesic. In the language of mechanics it measures the rate of relative acceleration 
of two particles moving forward on neighboring geodesies. 

The mixed scalar curvature is the following function, see [51 II 1 1 [IS]: 

Sc m i x (c/) = / / R(s a , ej, £ a , €.%), 

where {e^, £ a }i<n,a<p is a local orthonormal frame on TM adapted to D and Djr. If a distribution 
cither D or Djr is one-dimensional and a unit vector field N is tangent to it, then the mixed scalar 
curvature is simply the Ricci curvature Ric s (iV, N). On a foliated surface (M 2 ,g) this coincides 
with the gaussian curvature: Sc m i x (<?) = Kic g (N, N) = K(g). Recall the formula, see [IS] : 

Sc mix ( 5 ) = div(F + fl» + \H\ 2 + \\T\\ 2 - \\h || 2 + \H T \ 2 - HMI 2 . (2) 

Integrating ([2]) over a closed manifold and using the Divergence Theorem, we obtain the integral 
formula with the total Sc m i x (<7). Formula ([2]) yields decomposition criteria for foliated manifolds 
under the constraints on the sign of Sc m i x (g), see [15] and a survey in [8j. 

Example 1. For any n > 2 and p > 1 there exists a fibre bundle with a closed (n + p)-dimensional 
total space and compact p-dimensional totally geodesic fibers, having constant mixed scalar cur- 
vature. To show this, consider the Hopf fibration tt : S 3 — > S 2 of a unit sphere S 3 by great circles 
(closed geodesies). Let (F,g\) and {B,g2) be closed Riemannian manifolds with dimensions, re- 
spectively, p — 1 and n — 2. Let M = F x S 3 x B has the metric product g = g\ x 52 • Then 
7r : M — > S 2 x B is a fibration with a totally geodesic fiber F x S . Certainly, Sc m i x (<7) = 2. 

One may ask the question: Which foliations admit a metric with a given property o/Sc m j x (e.g., 
positive or negative)! 

3. Flows of metrics on foliations. We shall examine the question using evolution equations. 
A flow of metrics on a manifold is a solution gt of a differential equation dtg = S(g) , where the ge- 
ometric functional S(g) is a symmetric (0, 2)-tensor usually related to some kind of curvature. This 
corresponds to a dynamical system in the infinite-dimensional space of all appropriate geometric 
structures on a given manifold. Denote by A4(D, Djr) the space of smooth Riemannian metrics 
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on M such that Djr is orthogonal to D. Elements of A4(D, Dp) are called adapted metrics to the 
pair (D, Djr). The notion of the D-truncated (r, 2)-tensor field S (where r = 0, 1) will be helpful: 
S{X\,X2) = S(X±, X^)- The D-truncated metric tensor g is given by g(X\,X2) = g(X\,X2) 
and g(Y, •) = for all Xi € D, Y £ Dj. For flows of D-conformal adapted metrics we have 
S(g) = s(g) g, where s(g) is a smooth function on the space of metrics on M. 

Rovenski and Walczak [11] (see also [9]) studied flows of metrics that depend on the extrinsic 
geometry of codimension-one foliations, and posed the question: 

Given a geometric property (P), can one find an P -truncated flow dtg = S(g) on a foliation 
(M,P) such that the solution metrics gt (t > 0) converge to a metric for which P enjoys (P)? 

Rovenski and Wolak [12] studied D-conformal flows of metrics on a foliation in order to prescribe 
the mean curvature vector H of D. Recently, authors |10j studied D-conformal flows of metrics on 
a totally geodesic foliation. In aim to prescribe the sign of Sc m i x (g i ), we shall study the following 
mixed scalar curvature flow of metrics gt- 

d tg = -2(Sc m Ug)-<s>)g- (3) 

Here $ : M — > R is a leaf-wise constant function, its value is clarified in what follows. By Lemma[2l 
the flow ([3D preserves harmonic (in particular, totally geodesic) foliations. 

Example 2. (a) Let (M 2 ,go) be a surface of Gaussian curvature K, endowed with a unit geodesic 
vector field N . Taking $ = 0, we reduce © to the following form: 

d t g = -2K{g)g. (4) 

This PDE looks like the Ricci flow on surfaces, but its RHS includes the truncated metric g instead 
of g. Using the formula N(k) = k 2 + K for the geodesic curvature of orthogonal to N curves, 
([ID yields the Burgers equation dtk + N(k 2 ) = N(N(k)), see [TU]. For a surface of revolution 
M 2 C R 3 with standard parametrization r = [p(x)cos6>, p(x)sm9, h{x)\ (0 < x < I, \6\ < 7r), 
where p > and (p 1 ) 2 + (h') 2 = 1, the profile curves 9 = const are geodesies tangent to N. Since 
k = — iV(logp), we have dtp = N(N(p)), this provides geometrical interpretation of the classical 
relation (i.e., Cole-Hopf transformation) between the Burgers and the heat equations. 

(b) For the Hopf fibration it : (S 2m+1 ,g can ) — > CP m of a unit sphere with fiber S , the orthog- 
onal distribution D is non-integrable while h = 0. By Q, Sc m ; x (<7) = 2m. Thus, g ca _ n on S 2m+1 is 
a fixed point of the flow ([3D with $ = 2m. 

The basic question that we address in the paper is the following: Given a Riemannian manifold 
(M, g) with a harmonic foliation P, when do solution metrics gt of f3[) converge to the limit metric 
g such that Sc m i x (<7) is positive or negative! 

4. The nonlinear Schrodinger heat equation. The flow @ provides the Burgers equation 
for the mean curvature vector H of D, similarly to the flow Q in Example [2ja). If H is leaf-wise 
conservative, then the potential obeys the non-linear Schrodinger heat equation, see Proposition [2[ 
Under certain conditions the asymptotic dynamics of its solutions is the same as of the linear 
Schrodinger heat equation, when the corresponding Burgers type equation for H has a single-point 
global attractor. (The case of totally geodesic P with integrable D yields the forced Burgers 
equation for H, the asymptotic behavior of solutions and corresponding linear Schrodinger heat 
equation have been studied in [l0]). The flow admits a unique global solution g t , whose Sc m i x 
converges exponentially to a leaf-wise constant. The metrics are smooth on M when all leaves are 
compact and have finite holonomy group. Using rescaling along D one may provide convergence to 
metrics with Sc m i x as positive so negative. 

5. The structure of the paper. In Introduction we pose the question about prescribing the 
mixed scalar curvature of a foliation and define the flow of metrics depending on this kind of curva- 
ture. Sections [T] and [2] contain main results and their proofs. These use results of Section [3] about 
the non- linear Schrodinger heat equation on a closed manifold. Throughout the paper everything 
(manifolds, foliations, etc.) is assumed to be smooth (i.e., C°°-differentiable) and oriented. We 
also assume that all the leaves of a foliation P are compact minimal submanifolds. 
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1 Main results 



An important step in the study of evolutionary PDEs is to show short-time existence/uniqueness. 

Proposition 1. Let J 7 be a harmonic foliation on a closed Riemannian manifold (M, go)- Then 
the linearization of is the leaf-wise parabolic PDE, hence f3|) has a unique solution gt defined 
on a positive time interval [0, to) an d smooth on the leaves. 

We denote V^/ := (V/)^. Given a vector field X and a function (j) on M, define the functions 
using the leaf-wise derivatives: the divergence divj^X = Yla=i di^aX, e a ) and the Laplacian 
= divjr(V -7r (/>). Notice that V^, divj- and Ajr (i.e., along the leaves) are t-independent. 

Based on the "linear algebra" inequality n||/i|| 2 > \H\ 2 with the equality when D is totally 
umbilical, we introduce the following function (a measure of "non-umbilicity" of D): 

p D ■= n - 2 (n\\h\\ 2 - \H\ 2 ) > 0. (5) 

For p = 1, we have /3d = n~ 2 ^j<j(^ — %) 2 ' where k{ are the principal curvatures of D. 

The Schrodinger operator is central to all of quantum mechanics. By Proposition [3] (in Sec- 
tion [2?2]), the flow of metrics ([3]) preserves the leaf-wise Schrodinger operator TL, given by 

H{u) = -Ajru-(3 D u. (6) 

The spectrum of % is an infinite sequence of isolated real eigenvalues Ao < Ai < A2 < • • • < Aj < . . . 
(leaf- wise constant functions) counting their multiplicities, and lim^oo Xj = 00. One may fix in 
L2 an orthonormal basis of corresponding eigenfunctions {e^}, i.e., T-L(ej) = Xjej. 

Remark 1. If the leaf F(x) through x £ M is compact then Ai > Ao and the eigenfunction eo 
(called the ground state) for the least eigenvalue Ao may be chosen positive, see Proposition 01 

The so called fundamental gap Ai — Ao of T~L has interesting mathematical and physical impli- 
cations (e.g., in refinements of Poincare inequality and a priori estimates, it also can be used to 
control the rate of convergence of numerical methods of computation), see survey in pp. 

We shall say that a smooth function fit, x) on (0, 00) x F converges to f(x) as t — > 00 in C°°, if 
it converges in C fc -norm for any k > 0. It converges exponentially fast if there exists u > (called 
the exponential rate) such that lim^oo e ut \\f(t, •) — f\\c^ = for any k > 0. 

Define a function * := (\\T\\ 2 go + \\hjr\\ 2 gQ ) and the domain U = {x E M : ||T||g + ||Mlg > 0}. 

Proposition 2. Let a family of metrics gt (0 < t < to) solves (0) on a foliated manifold (M, J 7 ). 
Suppose that J 7 is a harmonic for go and Hq = — raV^log uq) for a function uq > 0. Then 

8 t H + V T g{H, H) = nV^div^tf) - n 2 V T (3 D + nV^GlTlg + \\hjr\\ 2 g ), (7) 

moreover, Ht = — nV^logu) for a positive function u : M x [0,to). 

(i) IfJJ = M then u = ^ 1 / 4 (||T||2 t + ||^||^ t ) _1/4 and the non-linear PDE is satisfied 

(l/n)d t u = A T u + (f3 D + $/n)u- (^/n)u- 3 , u(-,0)=u . (8) 

(ii) If U = then (J7J) reads as a forced Burgers equation 

d t H + V T g(H, H) = raV^(div T H) - n 2 V T (3 D (9) 
and the function u (i.e., the potential for H) may be chosen as a solution of the linear PDE 

(l/n)d t u = Ajru + l3 D u, u(-,0)=u . (10) 
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Under certain conditions ([9]) has a single-point global exponential attractor. Kirsch-Simon [6] 
studied the Burgers PDE on W 1 and proved the polynomial convergence of a solution. 
Based on Proposition [2Jii) and methods of [10], we obtain the following. 

Theorem 1. Let T be a totally geodesic compact foliation with integrable orthogonal distribution 
on a Riemannian manifold (M,go), and Hq = — nV" r (loguo) for a function uq > 0. Then |5p has 
a unique global solution gt (t > 0) smooth on the leaves. If $ = n \q then, as t — > oo, the metrics 
gt converge in C°° with the exponential rate n(Ai — Ao) to the limit metric g and 

Sc mix (<?) = nA , H = -nV- F (log e ). 

Moreover, if the leaves have finite holonomy group, then all gt and g are smooth on M. 

Example 3. a) Let (M, gt) = Mi x f t M2 be a family of twisted products, i.e. the manifold 
M = Mi x M2 with the metrics gt = ft9i + 92 where ft S C°°(Mi x M 2 ) are positive functions. 
The submanifolds {x} x M2 of a twisted product are totally geodesic, while Mi x {y} are totally 
umbilical with the mean curvature vector H t = — V^log/t), see [7]. Hence twisted products 
provide example with leaf- wise conservative Hq. The above metrics gt solve ([3]) with $ = if and 
only if the warping function obeys the heat equation ^ dtf = Ajf/. The conditions above yield the 
Burgers type equation for H (see ([9]) with (3d = 0) 

d t H + V T g(H,H) = nV F {A\w T H). (11) 

If Mi and M2 are closed then gt {t > 0) converge as t — > 00 in C°° with the exponential rate Ai to 
the metric g of the product (Mi,f 2 g 1 ) x (M 2 ,5 2 ), where f{x) =J M2 f(0,x,y)dy g , see [TO]. 

b) Let (M 2 ,go) be a two-dimensional torus with Gaussian curvature K, foliated by closed 
geodesies with a unit tangent vector field N. Suppose that the geodesic curvature k of orthogonal 
(to the leaves) curves obeys k = N(ipo) for a smooth function ?/>o on M 2 . The operator ([6]) coincides 
with the leaf-wise Laplacian — Ajr, hence Ao = and eo = const. We also have \& = T = (3d = 0. 
By Theorem [H the flow of metrics (UD on M 2 admits a unique solution gt (t > 0) converging as 
t — > 00 to a flat metric, and the leaves compose a rational linear foliation on the torus. 

The central result of the work is the following. 

Theorem 2. Let (M,go) be a Riemannian manifold endowed with a harmonic compact foliation J 7 . 
Suppose that U = M, Hq = — nV^logito) for a function uq > 0, and 

$ > -n(3 D , |nA + $| > max(||T||^ + \\hjr\\ 2 ) ( max(u /e ) / min(n /e )) 4 . (12) 

F F F 

Then |3p admits a unique global solution gt = 9t~ + 9 {t > 0) smooth on any leaf F, moreover, if 
all leaves have finite holonomy group, then gt are smooth on M . As t — > 00, we have 

Sc mix (3t) nA , H t -raV^log e ), HMU °> 

and convergence in C°° has the exponential rate na for any a € (0, min{Ai — Ao, 4 |Ao + $/n|}). 

Corollary 1. Let (M, 50) be a Riemannian manifold endowed with a totally geodesic compact 
foliation T . Suppose that U = M, Hq = — nV^logUo) for a function uq > and 

$>-n/3 D , |nA + ^| > max ||T||g (max(u /eo) /min(ii /eo)) 4 . (13) 

F F F 

Then the claim of Theorem @ holds. 

Let u obeys |8| ona leaf F. Define, see Section [3l uq = Uq + f °° qo(r) dr where Uq = j F uq eo dx 
and q (r) = _ e ( A o+*/^ f F (^/ n ) n" 3 e dx. 
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Corollary 2. Let (M,g) be a Riemannian manifold endowed with a harmonic compact foliation T . 
Suppose that U = M and H = — nV JF (loguo) for a function uq > 0. 

(i) If Xo < then there exists D-conformal to g metric g with Sc m i x (g) < 0. 

(ii) If Xo > — nCi^^Gl^llg + ll^-HIg) then there exists D-conformal to g metric g with Sc m i x (g) > 0. 

We consider applications to a Riemannian manifold (M,g) with a unit vector field iV (i.e., p = 1 
or n = 1). In this case, Sc m i x (g) is the Ricci curvature Kic g (N, N) in the iV-direction. 

Case p = 1. Let TV" be a unit tangent to a foliation by closed geodesies. Denote by t\ = Tr ff h 
the mean curvature of the orthogonal (to N) distribution D. We have \\hj?\\ = \\Hj?\\, and (j22[) 
reads as 

Ric g (N,N) = \\Tf g + N( n ) - (l/n) T ? -np D . (14) 

Suppose the condition n = — nN(\oguo) for a smooth function no > on M. Let the metric 
evolves as 

d t g = -2(m Cg (N,N)-<5>)g, (15) 
see (|3|), then t\ obeys the PDE along ^-curves, see ([7]), 

dtn + Ntf) =nN{N{T l ))-n 2 N{p D )+nN{\\Tf g ). (16) 

If D is integrable, then t\ = —n JV(log u) , where the function u( ■ , t) > may be chosen as a solution 
of the linear Schrodinger heat equation, see ()10p . 

(l/n)d t u = N(N(u)) + /3 D u, u(-,0)=u . (17) 

By Theorem [H the flow (|15|) admits a unique global solution gt (t > 0). We have convergence as 
t — > oo of gt — > g, T\ — > — n iV(logeo) and Ric St (Af N) — > nXo with the exponential rate n(Ai — Ao). 
If the foliation by iV-curves has finite holonomy group then the metrics gt and g are smooth on M. 

If D is nowhere integrable, then also t\ = —n N(logu), where the function u = 1 I/ 1//4 ||r|| 5t 1 ^ 2 
(with * := ul \\T\\ 2 go > 0) solves the non-linear Schrodinger heat equation, see ([8]), 

(l/n)d t u = N(N(u)) + (f3 D + $/n)u-(y/n)u- 3 , u(-,0)=u . (18) 

If conditions (fTHj) are satisfied, by Corollary [TJ the flow (fT5"j) admits a unique global solution 
9t (t>0). We have convergence as t — > oo of functions t\ — > — nA^logeo) and Ric St (iV, N) — > nXo 
with the exponential rate na for any a € (0, min{Ai — Ao, 4 |Ao + $/ra|}). If the foliation by 
A^-curves has finite holonomy group then the metrics gt are smooth on M. By Corollary [2l 

(i) If Ao < then there exists D-conformal to g metric g such that Ricg(A r , N) < 0. 

(ii) If Ao > — ^(f^) 4 ||^~1lg then there exists .D-conformal to g metric g such that Ricg(A r , N) > 0. 
Case n = 1. Let ./V be orthogonal to a compact harmonic foliation T of codimension one. Then 

T = /3 D = Hjr = 0, and fl22J) reads as 

Ric 9 (iV,iV) = diw T H-g{H,H) + \\h F \\ 2 g . (19) 

Suppose the condition Hq = — \7- F (log«o) for a smooth function uo > on M. 
If hjr ^ (general case) then = — V' F (logu) where, see dHJ, 

t u = A^u + (/3 D + $)u-^n" 3 , u(-,0)=«o. 

If the following conditions are satisfied: 

<I> > —nPo, \nXo + <&| > max ||/ij-||^ (max(uo/eo) / min(uo/eo)) 4 , 

.P i* 1 .F 

by Theorem [2j the flow (|15p admits a unique global solution <^ (t > 0). As i — > oo, we have 
convergence with the exponential rate a for any a € (0, min{Ai — Ao, 4 |Ao + 3>|}) 

F^-V^(logeo), Ric 5t (iV,;V) -> A , \\h T \\ gt ^0. 
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If the foliation T has finite holonomy group then the metrics gt are smooth on M. By Corollary [2 

(i) If Ao < then there exists L>-conformal to g metric g such that Ricg(./V, N) < 0. 

(ii) If Ao > — (i^) 4 ||^J"llg then there exists D-conformal to g metric g such that RiCg(iV, N) > 0. 
In case of totally geodesic foliation F, i.e., hj = 0, the operator ([6]) coincides with — Ajr (hence, 

Ao = and eo = const), and (|22|) reads as 

Ric 3 (iV, N) = div^ H - g(H, H). (20) 

Let the metric evolves by (|15p with $ = 0. By Proposition [21 H obeys the homogeneous Burgers 
equation (jlip with n = 1. Suppose that the curvature vector H of iV-curves is leaf- wise conservative. 
Substituting H = — W^ilogu) we find that the function u > obeys the heat equation dtu = Aj-it. 
The global solution of above PDE on the leaves obeys u := lim u(t,x) = j F uq(x) dx/Vol(-F x , g). 

Since eo = const, we have H = lim H(t, •) = 0. Then Ricg(iV, N) = 0, where g = lim gt. 

t— >oo t— >oo 

2 Proof of main results 
2.1 Preliminaries 

The Levi-Civita connection V of a metric g on M is given by the formula 

2g(V x Y,Z) = X(g(Y,Z)) + Y(g(X,Z)) - Z(g(X,Y)) 

+ g([X,Y],Z)-g([X,Z},Y)-g([Y,Z},X) (X,Y,ZeTM). (21) 

The co-nullity operator is defined by Cjy(X) = — (VxA r )- L , for X £ D, N G Z?jr- One may 
decompose C into symmetric and skew-symmetric parts as Cn = An + T^. The Weingarten 
operator An of D and the operator are related with tensors h and T, see ([T]), by 

g(A N (X),Y)=g(h(X,Y), N), g(T N (X), Y) = g(T(X,Y), N), X,Y € D. 

Next lemma represents Sc m i x (g) in a more convenient than ([2]) form. It shows that ([3]) can be 
reduced to the leaf- wise evolution equation (with the space derivatives along T only). 

Lemma 1 (see also |10j). Let T be a harmonic foliation on (M,g). Then 

Sc mix (g) = \\T\\ 2 g + UMIg + div T H- g(H, H)/n -nfo. (22) 

The notion of holonomy uses that of a germ of a locally defined diffeomorphism (i.e., an equiva- 
lence class of certain maps). The germs of diffeomorphisms (M n , 0) — > (M n , 0) with fixed origin form 
a group, denoted by Diffo(M n ). We denote by Diff(j~(]R n ) the subgroup of germs of diffeomorphisms 
which preserve orientation of M n . Let (M, J-) be a foliated manifold, x,y be two points on a leaf 
F and Q x ,Qy be transversal sections (diffeomorphic to K n ). To any path a from x to y in F we 
associate a germ of a diffeomorphism h(a) : (Q x ,x) — > (Q y ,y) called the holonomy of the path, 
for x = y, this is a generalization of a first return map. We obtain a homomorphism of groups 
h : 7ri(F) — > Diff (M n ). The image Hol(F) := h(-K\{F)) is called the holonomy group of F. If T is 
transversally orientable then Hol(-F) is a subgroup of Diff^M"). 

Certainly, Hol(F) is finite when the first fundamental group, ttx(F), is finite. Note that a foli- 
ation whose leaves are the fibers of a fibre bundle has trivial holonomy group. 

Definition 1. A foliated bundle is a fibre bundle p : E — > F admitting a foliation T whose leaves 
meet transversely all the fibers E x = p~ 1 (x) (x € F), and the bundle projection restricted to each 
leaf F' € T is a covering map p : F' ^ F. There is a representation h x : 7Ti(i ? , x) — > Diff(£' x ) 
of the first fundamental group tti(F, x) in the group of diffeomorphims of a fiber, called the total 
holonomy homomorphism for the foliated bundle. 
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Local Reeb stability Theorem (see [31 Vol. I]). Let F be a compact leaf of a foliated manifold 
(M,J-) and Hol(-F) is finite. Then there is a normal neighborhood pr : V — >• F of F in M such that 
(V, J-\y,pr) is a foliated bundle with all leaves compact (and transversal to the fibers). Furthermore, 
each leaf F' C V has finite holonomy group of order at most the order o/Hol(-F) and the covering 
pri^/ : F' — > F has k sheets, where k < order of Hol(-F). 

In other words, for a compact leaf F with finite holonomy group, there exists a saturated 
neighborhood V of F in M and a diffeomorphism from E := F x Hol ^ W 1 under which T\y 
corresponds to the bundle foliation on E. Here F is a covering space of F associated with Ho\(F). 

The following method construction foliations is related to local Reeb stability Theorem. Suppose 
that a group G acts freely and properly discontinuously on a connected manifold F such that 
F/G = F. Suppose also that G acts on a manifold Q. Now form the quotient space E := F xq Q, 
obtained from the product space E := F x Q by identifying (gy, z) with (y, gz) for any y £ F , g € G 
and z 6 Q. Thus E is the orbit space of E w.r.t. a properly discontinuous action of G. It is also 
Hausdorff, so it is a manifold. The projection pr x : E — > F induces a submersion pr : E — > F, so 
we have the commutative diagram (F — >• F) o pr x = pr o {E — » S). The map pr has the structure 
of a fibre bundle over F with vertical fiber Q. (Fibre bundles which can be obtained in this way 
are exactly those with discrete structure group). We claim that E admits also horizontal leaves, 
so that pr maps each leaf to F as a covering projection. Indeed, the foliation T on E, which is 
given by the submersion pr 2 : E —¥ Q, is invariant under the action of G, and hence we obtain the 
quotient foliation T = F/G on E. If z € Q and G z C G is the isotropy group at z of the action 
by G on Q, then the leaf of E obtained from the leaf F x {z} is naturally diffeomorphic to F/G z , 
and pr restricted to this leaf is the covering F/G z — > F. 

2.2 D-conformal adapted variations of metrics 

Let gt be a smooth family of metrics on N and S = dtg. Since the difference of two connections is 
a tensor, c^V* is a (l,2)-tensor on (M,g t ). Differentiating (|2T]) with respect to t yields, see 

2g t ({d t V t )(X,Y),Z) = (V$rS)(y,Z) + (V l Y S)(X,Z) - (V%S){X,Y) (23) 

for all t-independent vector fields X, Y, Z € T(TM). If S = s(g)g, for short we write 

d t g = sg (24) 

for a certain t-dependent function s on M. In this case, the volume form dvol evolves as [3] 

(d/dt)(dvol t ) = (n/2) s t dvol t . (25) 

Lemma 2. For D-conformal adapted variations (|24p o/ metrics we have 

dthjr = -a /or, d t H T = -sH r . (26) 

Hence, variations (|24p o/ metrics preserve harmonic and totally geodesic foliations. 

Proof. Let (i > 0) be a family of metrics on (M, J 7 ) such that dtgj = #(<?), where the tensor 
5(y) is D-truncated. Using (|23|) , we find for X £ D and £,77 6 Dj-, 

2ffi(«,'))J) = ^(Vl^ + c^V^), ^) = (V|5)(X,r ? ) + (V*5)(X,0-(V t x 5)(e,r?) 
= -5(V|r?,X) -5(V*^,X) = -2 5(M^»/)^)- 

Assuming S"(g) = s(g)g, we have f|26[) 1 . Tracing this we have ()26p 9. By the theory of ODEs, if 
Hjr = or hjr = at t = then respectively Hjr = or hjr = for all t > 0. □ 

The proof of the next lemma is based on ([23]) with S = s g. 
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Lemma 3 (see [TO] and |12j). For D-conformal adapted variations (|24|) of metrics we have 

d t A N = ~N(s)id, d t Tl = -sT i N d t H = ~V T s. (27) 

where N E Dj. Hence, variations ()24|) of metrics preserve conformal foliations. 

Note that variations (|24|) of metrics preserve the domain U. Denote 5t := 11^11^ + Harlot > 0. 

Proposition 3 (Conservation laws). Let gt (t > 0) be D-conformal adapted metrics (|24|) on 
a foliated manifold (M,J-,D) such that Hq = — nV"^(loguo) for a positive function uq 6 C°°(M). 
Then the function (3d and the vector field Ht — ?V (log<5t) on U don't depend on t. 

Proof. Using Lemma and g(H, • ) = 0, we calculate 

dt\\h\\ 2 = Tr(,4 2 J = 2^ Tr(A ea d t A ea ) = -J2 e tt (s) Tr A £a = -g(Vs, H), 

d t g(H,H) = S g(H,H) + 2g(d t H,H) = -ng(Vs,H). 

Hence, ndt^D = dt \\h \\ 2 g — ^ dtg(H,H) = 0, that is the function /3d doesn't depend on t. 
For any function / 6 C (M) and a vector iV € Djr, using (dtg)(-,N) = 0, we find 

g(V^(d t f),N) = N(d t f) = 8 t N(f) = d t g(y T f,N) = g(d t (V T f), N). 

Hence V T {d t f) = d t {V T f). By LemmaEl we find 

9t ||T || 2 = -5, ^ Tr ((Tjj 2 ) = -2 ^ Tr (3^^) = 2 S ^ Tr ((Tjj 2 ) = -2 S ||T || 2 . 

(28) 

Similarly, using the proof of Lemma [5J we obtain 

d t \\h T \\ 2 = -2s\\hrf. (29) 

By (f28|) and (f29|) . the domain J7 is t- independent. On U we have dt(logdt) = —2s. (If T = or 
/ijr = at x € M then respectively ||r(x)|| 9t = or ||/i^r(x)|| 9t = for all t > 0). For x £ U, using 
V^dt = dtV^ , we obtain 

Sfcflt = (n/4)d t V^(log5 t ). 
Hence the vector field — (n/4)V" F (log <5f) (defined on J7) doesn't depend on t. □ 

2.3 Proofs of Propositions [TH2], Theorem [2] and Corollary [2] 

Proof, (of Proposition [1]) Let gt (t > 0) be a family of D-conformal metrics on a foliated manifold 
(M,J-), i.e., dtgt = stgt- By Lemma[21 T is harmonic with respect to all gt- We differentiate (|22|) 
by t, and apply d t V T H = V T d t H, ([25]) . (f29|) and Lemma [3] to obtain 

d t Sc mix (<7 t ) = -(n/2) Ajr s + ff (Vs, iJ) - 2s 5*. (30) 

Next, consider D-conformal metrics g T = go + rs 50 f° r small r£l, where s = — 2(Sc m i x (g , o) — 
By ()30p . the linearization of (|3j) at go is the linear PDE on the leaves 

<9 T Sc mix (£ r )| r= o = nAjrSc^igr) -2g (VSc mix (g T ),H(g )) 
- 2 (Sc mix (g ) - $o - $)(Sc mix (s T ) - $) 

and Sc m i x (go) is bounded. The result follows from the theory of linear parabolic PDEs (see [2]) 
and the finite holonomy assumption (i.e., the local Reeb stability Theorem in Section [2.1|) . □ 
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Proof, (of Proposition [2]) By Proposition[TJ ([3]) admits a unit local smooth solution gt (0 < t < to). 
The functions Sc m ; x (^), H t , \\T\\ gt and ||/ij-|| St etc. are then uniquely determined for < t < to- 
From d27|)3 with s = -2 (Sc mix (#) - $) and using (J22J) we obtain 0. 

(i) By Proposition [3jii) , Ht — (n/4)V' ?r (log 5t) = X for some vector field X on M. Since iio is 
conservative, X = — (n/4)V' 7r (log for some smooth function $ > on M. Hence, 

H t = -nV T log (^/% 1/4 ) 
and ^ := Uq 5o by conditions. Define a smooth function u := VP 1 / 4 ^ 1//4 onf/x [0, to) an d calculate 

d t {\og5 t ) = -Ad t \og(8; 1/A ) = -4d t log(*- 1/4 u) = -Ad t {\ogu). 
From the above and (j22|) we then obtain 

^(logtt) = -(l/A)d t (\og5 t ) = -Sc mix (9t) + $ 

= nA J r(logu) + nfl'(V- 7r (logn), V jr (logu)) +n/3 D + $ - ^n~ 4 . 

Substituting 

<%(logu) = u^dtu, V^logu) = it -1 V- 7r «, Aj-(logu) = u -1 Ajr it - vT^g^u, V^it), 
we find that u solves the non-linear Schrodinger heat equation ([8|). (ii) See the proof in [10]. □ 

Proof, (of Theorem[2]) (i) By Proposition [IJ there exists a unique local solution gt on M x [0,io)> 
and by Proposition [2)^ ii) there is corresponding positive function u satisfying ([8]) with it(-,0) = uo 
such that H = — nV" F (logu) and (J7J) holds. Note that conditions (|12p yield 

|nA + $| > max{(5o(uo/eo) 4 }/min(uo/e ) 4 <^ (u^) 4 > {^ + /n)/\\ + $/n| 

F F 

that is equivalent to ([55]) with /3 = /fo + $/n > 0. The quantities Uq and ^ + are defined in ([55]) . 

By Theorem [3j one may leaf-wise smoothly extend a solution of (JSJ) on M x [0, oo), hence Ht(x) 
is defined for all t > and is smooth on the leaves. 

By TheoremSfi), u(- ,t) — > oo as i -4 oo with the exponential rate na for any a £ (0, min{Ai — 
Ao, 4 | Ao + 3>/n|}). Hence, the function <5j = ^u -4 is leaf-wise smooth, moreover, ||T|| gt — > and 
II^HIst — >■ as i — )■ oo. 

By Theorem Hl^ii) , Ht = — nV~ F (logit) approaches in C°°, as t — > oo, to the vector field H = 
— nV"^(logeo) and div Ht approaches to the smooth function — nAjr(log e o)- Recall that 

-A T eo - (Pd + $/n)e = (A + $/n)e . 

Hence, divff t — <7(-^t) H t )/n approaches to — n(Aj-eo)/eo = nAo + n /?d as t — > oo. Using ([22]) . we 
find a smooth on the leaves function Sc m i x (x,i) which approaches exponentially to nAo as t — > oo. 
The leaf-wise smooth solution to (J3j is given by g t = go exp(— 2 J (Sc m i x (x, r) — dr) for t > 0. 

(ii) The smoothness of (ft on M follows from the finite holonomy assumption and results of 
Section 13.41 (see also |10j). Indeed, let F be a leaf. By the local Reeb stability Theorem (see 
Section [2. ip . there is a normal neighborhood pr : V — > F of F (with a smooth normal section Q - 
an open n-dimensional disk) such that (V, J-iy,pr) is a foliated bundle diffeomorphic to F Xn \(F)Q- 
There is a regular covering pr : F — > F with covering group Hol(-F), since this group is finite, F is a 
compact manifold. The normal neighborhood V, as a fiber bundle over F, can be pulled back via pr 
to a bundle V over F with the same fiber. The standard pull-back construction yields a canonical 
covering pr : V — > V (a submersion), enabling us to lift the foliation F\y to a foliation T of V, 
transverse to the fibers and having F as a leaf. Let ieF with pr(x) = x. Since the covering group 
is exactly the holonomy group of F, pr* : m(F,x) — > 7ri(.F, x) injects n±(F, x) onto the subgroup 
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of tti(F, x), and the leaf F of J- has trivial holonomy. By [31 Theorem 2.4.1], F has a neighborhood 
in V that is a foliated disk bundle with all leaves diffeomorphic to F. Let (ft be the lifts of the 
metrics gt from V onto the product V . The corresponding foliation on V is harmonic, the lift of 
the laplacian and potential function, /3d, smoothly depend on q on V. (In case of totally geodesic 
foliation the leaves are isometric one to another, see for example [5]. Hence the leaf- wise laplacian 
on V doesn't depend on q £ Q). The functions e"o and solution to (j5|) . Uf (t > 0). are smooth 
on V" (where Q can be replaced by smaller normal section), and they are lifts of leaf- wise smooth 
functions eo and Ut (t > 0) on V or on a smaller neighborhood W of F. The vector fields Ht and 
H are smooth on V, and they are lifts of smooth vector fields Ht and H on a neighborhood of .F. 
By the above, gt are smooth on M. □ 

Proof, (of Corollary [2]) The metrics of Theorem [2] diverge as t —¥ oo with the exponential rate 
H = n\ + 

3C> 1, VX G D, Vi > : C^e'^g^X^X) < g t (X,X) < Ce~^g (X,X). 

Consider D-conformal metrics gt = e~ fJ-t g^- + g. By (|27p ^. Ht = Ht. The function v = e flt u 
converges as t — > oo to uoeo, where «o = (u,eo)o = + Jo°° 9o( r ) dr, see Theorem H]and l[67|) . 

Denote 6 t := ||T||| 4 + ||^||? . We then have v 4 (-,t)5 t = * for all t. Hence, S t = e~ 2 ^'^ 
converges as t — > oo to a smooth function «q/(«oCo) 5o> and the metrics gt converge as t — > oo to 
the metric g^ = (u e /u ) 2 gQ + g. By ([22j), Sc mix (^) = Sc mix (^) -6 t + S t converges as t -)• oo to 

Sc mix (goo) = nA + «q/ (n e ) 4 5o- □ 

3 The non-linear Schrodinger heat equation 

The section plays an important role in proofs of main results (see Section [2]). 

Let (F, g) be a closed p-dimensional Riemannian manifold, e.g., a leaf of a foliation T . Functional 
spaces over F will be denoted without writing (F), for example, L2 instead of L2(F). 

Let H l be the Hilbert space of differentiable by Sobolev real functions on F, with the inner 
product ( •, • )i and the norm In particular, H° = L2 with the product ( •, • )o and the norm || - ||o- 

If E is a Banach space, we denote by || • \\e the norm of vectors in this space. If B and C 
are real Banach spaces, we denote by B r (B,C) the Banach space of all bounded r-linear operators 
A ■ Wi=i B -> C with the norm \\A\\ B r {B ^ c) = sup Vu _ yVrGB \ ■ If r = 1, we shall write 

B(B, C) and A(-), and if B = C we shall write B r (B) and B(B), respectively. 

If M is a /c-regular manifold or an open neighborhood of the origin in a real Banach space, 
and N is a real Banach space, we denote by C k (M,N) (k > 1) the Banach space of all C fc -regular 
functions / : M — > N, for which the following norm is finite: 

\\f\\c k (M,N) = sup a;eM max{||/(a:)||^, max!< r < fc \\d r f(x)\\ B r {TxMtN) }. 

Denote by || • \\ck, where < k < 00, the norm in the Banach space C k ; certainly, || • \\c when 
k = 0. In coordinates (x±, . . . ,x p ) on F, we have ||/||cfc — m&x-xeF max| a |< fc \d a f(x)\, where a > 
is the multi-index of order \a\ = Y2i=i a * an d d a is the partial derivative. 

3.1 The linear Schrodinger heat equation 

For a smooth (non-constant in general) function /3 on (F,g), the Schrodinger operator, see ([6]), 

n(u) := -Au- f3u (31) 

is self-adjoint and bounded from below (but it is unbounded). The domain of definition of % is H 2 . 
The spectrum o~(7i) of H (with fixed C) consists of an infinite sequence of isolated real eigenvalues 
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Ao < Ai < . . . Xj < ... counting their multiplicities, and Xj — > oo as j — > oo. If we fix in L2 
an orthonormal basis of corresponding eigenfunctions {ej} (i.e., H(ej) = Xjej) then any function 
u E L2 is expanded into the series (converging to u in the Z^-norm) 

u(x) = . ^ Cj ej(x), Cj = (u,ej)o = J u(x)ej(x)dx. (32) 

The proof of (|32p is based on the following facts. Since by the Elliptic regularity Theorem, we have 
Ti^ 1 : L2 — > H 2 when and the embedding of H 2 into L2 is continuous and compact, see [2], then 
the operator %~ l : L2 — > L2 is compact. This means that the spectrum <r(T-L) is discrete, hence by 
the spectral expansion theorem for compact self-adjoint operators, {&j}j>o form an orthonormal 
basis in L2, see [5j I, Ch. VII, Sect. 4; and II, Ch. XII, Sect. 2]. One can add a constant to j3 such 
that 71 becomes invertible in L2 (e.g., Ao > 0) and Ti" 1 is bounded in L2. 

Proposition 4 (see jlUj). Let (3 be a smooth function on a closed Riemannian manifold (F,g). 
Then the eigenspace of operator (|31 j) . corresponding to the least eigenvalue, Xq, is one- dimensional, 
and it contains a positive smooth eigenf unction, eo- 

The following facts will be used. 
Sobolev embedding Theorem (see [2]). If a nonnegative k € 7L and I £ N are such that 
21 > p + 2 k, then H l is continuously embedded into C k . 

Elliptic regularity Theorem (see [2]). // H is given by (j3Tj) and ^ o-(H), then for any 
nonnegative k € Z we have W^ 1 : H k — > H k+2 . 

The Cauchy's problem for the heat equation with a linear reaction term has a form 

dtu = —T-L(u), u(-,0)=no, (33) 

We shall use the following scalar maximum principle [H Theorem 4.4]. 

Proposition 5. Suppose that X{t) is a smooth family of vector field on a closed Riemannian 
manifold (F,g), and f G C°°(F). Let u : F x [0, T) -4 R be a C°° supersolution to 

d t u > A g u + (X(t), Vu) + f(u, t). 

Suppose that ip : [0, T] — > R solves the Cauchy's problem for ODEs (d/dt)<p = f((j)(t),t), <^i(0) = C . 
If u(-,0) > C then u(-,t) > <p(t) for all t € [0, T\. (Claim also holds with the sense of all three 
inequalities reversed). 

3.2 Long-time solution 

Consider on (F,g) the Cauchy's problem for the non-linear Schrddinger heat equation: 

d t u = -H(u) - V(x) u -3 , u{x, 0) = u (x), (34) 

where T~i is defined in (|3ip . fi(x) > 0, and fy(x) > for any x 6 F. Note that after scaling the time 
(by n) and replacement VP/n — > * and [3d + $/n P, © reduces to ([34"]) . By [2j Theorem 4.51], 
the parabolic PDE ()34|) has a unique smooth solution u(-,t) for i € [0,to) an d some io > 0. 
Let f3(x) > j3~ > 0, hence (/3(x)u, u)q > [3~(u, u)q. Thus, for any u € Dom(%) we have 



(Hu, u) = J (\Vu(x)\ 2 - /3(x)\u(x)\ 2 ) dx < J (\Vu(x)\ 2 - P~\u(x)\ 2 ) dx = (-Au - /3~u, u) . 

Since (3~ is the maximal eigenvalue of the linear operator A + f}~ id (id is the identity function), 
by the variational principle for eigenvalues, we obtain that Ao < — /3~ < 0. 
Define the quantities 

W + = max^— , W =mm^— — , = max — — , u Q = mm — — . (35) 

x&F el[x) xeF eQ{x) x&F eo{x) xeF eo(x) 
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Lemma 4. Let u{x,t) > be a solution in the cylinder Ct = F x (0,T) of the Cauchy's problem 
\34\ ) with the condition, see (|12p . 

K") 4 > tt+/|Ao|. (36) 
Then the following a priori estimates are valid: 

w-(t) < u(x,t)/e (x) < w + (t) ((x,t) G C T ), (37) 

where 

W-(t)=(((u^-* + /\\ \)e- iX ° t +^/\\ \)K w + (t)=(((u+)^-/\X \)e-^ t +^-/\X \)\ (38) 
Proof. Since eo(x) > on F, we can change the unknown function in (j34[) : 

u(x,t) = eo(x)w(x,t). (39) 

Substituting into (|34p . we find 

e d t w = e Aw + (Ae + (3(x)e )w + 2 g(Ve , Vw) - e^ 3 (x)^(x)w~ 3 . 

Since Aeo + /3(x)eo = — Aoeo, we obtain the following Cauchy's problem for w(x,t): 

d t w = Aw - X w + 2#(V(loge ), Vw) - eQ 4 (x)^(x)w~~ 3 , w(-, 0) = u /e . (40) 

Let w(x,t) be a positive solution of the Cauchy's problem (|40p in a cylinder Ct- Then, using 
(|55b.b) . we obtain the differential inequalities for (x,t) G Ct- 

Aw - X w + 2 g(V(log e ), Vw) - ^ + w~ 3 < d t w < Aw - X w + 2 g(V(log e ), Vw) - ^"tt;" 3 . 

Using the scalar maximum principle of Proposition [5] and (|35b .d). we conclude that 

W-(t) < w(x,t) < w + (t) for any (x,t) G Ct, 

where W-(t) and w+(t) are the solutions of the following Cauchy's problems for ODEs: 

(d/dt)w- = -XqW- - $> + wZ 3 , W-(0) = Uq, 
(d/dt)w + = -X w + - $>~w+ 3 , w + (0)=Uq. 

One may check that these solutions are expressed by ([38]) , Note that if \E' = then (f3T|) reads as 
e~ Xot < u(x,t)/e (x) < uJe~ A °*. □ 



It will be convenient for us to change the unknown function in ([31 

v(x,t) = e Xat u(x,t). (41) 
Then we obtain the following Cauchy's problem for v(x,t): 

8 t v = -(H- X id)v -f?)- 3 e 4Ao ', v(x,0) = u (x). (42) 
Using dHJ, (|3"7|) and ([38]) . we estimate the solution v(x,t) of this problem when (f36|) holds: 

v- < v(x,t)/e (x) < v + ((x,t) G F x [0,oo)), (43) 
where the constants are given by 

^ = ((n -) 4 -^ + /|A |)^ v + =u+, (44) 
13 



Lemma 5. Let /3(x) > 0,^(x) > 0, u (x) > belong to C°°(F) and (J36D is satisfied. Then 

(i) for the partial derivatives by x of a positive solution v(x,t) of \J^%fy the estimates are valid: 
for any multi-index a = (a±, c*2, ■ ■ ■ , a p ) there exists C a > such that 

\d£v(x,t)\ < C a (l for any (x,t) G = F x [0,oo); 

(ii) for the partial derivatives by x of the function Q(x,t) = — ^(x)v ~^(x, i)e 4A °* we have the 
estimates: for any multi-index a there exists Q a > such that 

\d%Q{x,t)\ < Q a (l + t)\ a \e 4Xot for any (x,t) eC^. (45) 

Proof, (i) Denote for the sake of brevity Dj = d Xj (j = 1,2, ... ,p), D a = <9" = D ai D a2 ■ ■ ■ D ap . 
Differentiating both sides of (|42|) by the variables x\,... ,x p , we obtain the following equations for 
the functions p a (x,t) := d^v(x,t): 

d tPj = -(n - A id ) Pj + a{x, t) Pj + bj(x, t) + Cj {x, i)e 4A °' (j = 1, 2, . . . ,p), (46) 
dtPjk = —(H - A id )p jk + a(x, t)p jk + b jk {x, t) + c jk (x, t)e 4A °* {j, k = l,2,...,p) (47) 



and so on, where 



a( x ,t) = 3^e 4Aoi , bj(x,t) = Dj(/3)v, (48) 



v 



Cj (x,t) = -DjW v~ 3 , b jk (x,t) = Djiftpk + D jk {P)v, (49) 
c jk (x,t) = SDkWv^pj + SDjWv^Pk-D^v- 3 -12*v- 5 Pj p k (50) 

and so on. Let us make the change of unknown functions in equations (|46|) . (|47|) and so on: 

pj = pj e , Pjk = Pjk eo (51) 

and so on. Then in the same manner, as (|40p have been obtained from ()34p . we get 

d t pj = Apj + 2 <?(V(log e ), Vpj) + apj + bj/e + ( Cj /e ) e 4Aoi (j = 1, 2, . . . , p), (52) 
dtpjk = Afijk + 25-(V(log e ),Vp jk ) + ap jk + b jk /e + (c jk /e ) e 4A °* (j, fe = 1, 2, . . . ,p) (53) 

and so on. From (j4~8j) . (j4*3j) and (f5Tj) - (j5"2j) we get the differential inequalities 

Apj + 2 5 (V(log e ), Vpj) - a + \pj\ - bf - 6+e 4A °* < 

< ^ < Apj + 2 5 (V(log e ), Vjij) + a+ \p | + bf + c+e 4A °* 

for j = 1, 2, . . . ,p, where 

a + = 3^ + fZ 4 e 4Ao *, bf = pfv+, cf = ^fvZ 3 , 0+ = max \Dj(3\, $f = maxd-D^eo" 4 ). 

(54) 

Then by the maximum principle of Proposition [5] the estimate is valid: \pj(x,t)\ < pf{t) for any 
(x,t) E Cqo, where pf{t) is the solution of the Cauchy's problem for the ODE: 

(d/dt)pf = a+(t)\pf\ + bf + 6+e 4Aoi , p+(0) = p° (with p°j = max F 0)| ). 
As is known, 

p +(t) = f. exp ( f a+(r) dr) + T exp ( f a+(r) dr) (6+ + c+e 4A ° s ) ds 

JO JO is 

Since in view of (f54"|) i 

;>oo 

/ a+(r)dr<oo, (55) 
Jo 
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these circumstances mean that for any j € {1, 2, . . . ,p} there exists Cj > such that 

\Pj(x,t)\<Cj(l + t) ((M)eCoo), (56) 

hence, in view of (|51|) . we obtain claim (i) for \a\ = 1. 

Similarly we obtain from ([53"]) . ([5TJ]) . ([ST]) . ([35]) and ([51j]) the differential inequalities 



Ap ifc + 2 fl (V(loge ), Vp jfc ) - a + (t)|p jfe | - 6+ (i) - c+ (t)e 4A °* < d*^ 

where 



< Ap ifc + 2 5 (V(log eo), Vpjk) + a + (t)\p jk \ + b+(t) + c+(t)e 4Xot {j, k = 1, 2, . . . ,p), 



b %^) = +t) + /3+ k v+, P jk = max f \D jk (3\, *+ = max F | Z? jfc ^ | 4 , (57) 

= 3 (^t v - 4 + C'^>- 4 ) (1 + *) + 12 CA* + tC 5 (l + tf + «Z 3 . (58) 

Using again the maximum principle, we get the estimate for any (x,t) 6 C^: 

|PiJfe(x, t)| < P° fc exp ( jT* a+(r) dr) + jT* exp ( ^ o+(r) dr) (&+(*) + c+ (s)e 4A ° s ) ds, 

where f- k = max F |p jfc (-, t)\. Then d5T5J), (J5ZD and the fact that by ([SS]), J °° 6+ fc (s)e 4A ° s ds < oo 
imply that for any j, € {1, 2, . . . ,p} there exists Cj k > such that 

\Pjk{x,t)\ < C jk (l+t) 2 , (x,t) GCoo. 

Hence, in view of (|51j) . we obtain claim (i) for \a\ = 2. By induction we can prove (i) for any \a\. 
(ii) Estimate (gSJ) for |a| = 0, 

\Q(x, t)\ < (max F ^)tC 3 e 4A °*, (59) 

follows immediately from fj43j) . Estimate ([45 p for |a| = 1 and |a| =2 follows from claim (i), 
estimates (|4"5[) and equalities 



Dj(Q) = (3^- 4 ^ i ( W )- J D,(^)^ 3 )e 4Aot , 

^ifc(Q) = (3*«- 4 £> ifc (v) - D jk (*)v- 3 +3D k (*)v- 4 D j (v) +3D j (*)v- 4 D k (v) 
- 12 ^ v~ 5 Dj(v) D k {v))e 4Xot . 

For any \a\ claim (ii) can be proved by induction. □ 

Theorem 3. The Cauchy's problem {3$ with the condition (|36p admits a unique smooth solution 
u(x,t) > in the cylinder Coo — F x (0, oo). 

Proof. The positive solution u(x, t) of the Cauchy's problem (|34|) satisfies a priory estimates (|37p 
on any cylinder C^ where it exists. Then by the standard arguments, using the local theorem of 
the existence and uniqueness for semi-flows, we obtain that this solution can be uniquely prolonged 
on the whole cylinder Coo- Then, by Lemma (5[ all partial derivatives by x of u(x,t) exist in Coo- 
Hence, u is smooth on Coo- □ 
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3.3 Asymptotic behavior of solutions 

Theorem 4. Let u > be a smooth solution on Coo of the Cauchy's problem \S4-\j , see Theorem^ 
and (|36p is satisfied. Then there exists a positive solution u on Coo of the linear PDE 

d t u = An + (/3(x) + X )u (60) 

such that for any a S (0, min{Ai — Ao, 4 |Ao|}) 

(i) u = e~ Xot (u + 9(x,t)), where \\0(-,t)\\ C k = 0(e~ at ) as t -> oo for any k € N; 

(ii) V(logn) = V(logeo) + 6\{x,t), where ||0i(-, t)\\ C k = 0(e~ at ) as t — >■ oo for any k G N. 

Proof, (i) Let Gn(i, x, y) be the fundamental solution of (|60p (called the heat kernel). As is known, 

G (t,x,y) = J2.e-i* t e j (x)e j (y), (61) 

where (ij := Xj — Ao- Due to the Duhamel's principle, the solution v = e X(>t u, see (jSJ, of the 
Cauchy's problem (|42p satisfies the nonlinear integral equation 



v(x,t) = J^G {t,x,y)u (y)dy + ( jf G (t - r, x, y)Q(y, r) dy ) dr, (62) 

where Q is defined in Lemma [5] One may expand the functions v, uq and Q into Fourier series by 
the eigenfunctions ey. 

Vj(t)ej(x), Vj(t) = (v{-,t),ej) = v(y,t)ej(y) dy, (63) 
^o(^) = X^o^ 6 ^' u °j = ( u o,ej)o = J u (y)ej(y)dy, (64) 
QOM) = ^2. =Q <lj{t) e j( x ), Qj(t) = (Q(-,t),ej)o = / Q(y,t)ej(y) dy . (65) 
Then we obtain from (|6ip . ([62]): 

«,•(*) = uje - "** + ^ e-H^qjir) dr (j = 0, 1, . . . ). (66) 
In view of ([4"3"]) and ([55]) . we have 

/ |go(r)|dr< / (Q(., r ),eo)odr<— / e~ 4/3r dr < oo. 
Jo Jo ui Jo 

Substituting Vj(t) of ([66]) into ([63]) . we represent i> in the form v = u + 9, where 

oo r°° 

+ J2 -_i «?e~ w S. «o = «o + / ?o (r) dr, (67) 
J- 1 Jo 

/■oo \ oo 

jf gd(r) drjeo + Vj{t)ej, Vj(t) = J (r) dr 

Recall that uo = 0. Observe that u is a solution of ()60p with the initial condition 

u(-,Q)=u +(J g (r)dr)e . 
Let us estimate the terms of ([67]) and ([68]) . By ([4*3]) and ([35]) . we have 

/oo /-oo ,t,+ poo ^r,+ 

go(r)dr <y (|Q(., r )|,e )odr<-3-y e 4 ^ dr = ___ e 4A ^ 



u = ti e 



(68) 



4 An kv 
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Let us take k E N, I = + k/2\ + 1 and 7 < Ao- Using assumption uq € C°°(F) and the fact 
that Q(-,t) £ C°°(F) for any i > 0, we may consider the functions 

t«o := (ft - 7 id )W ^0, *) == (W - 7 id ) Z Q(-,i), 
which have the same properties: wq G C oc (F) and P(-,t) S C°°(F) for any f > 0. Let us represent: 

(it ,ej) ej = ((ft - 7id)" / ^ ,e J ) ej = (w , (ft - 7 id )" / e J ) 

= (wo,ej)o €j _ = (^-7id)~ / ((w ,e j )oe J ). 
Aj 7 

Similarly, we obtain (Q(-, t), ej)o &j = (ft — 7 id ) _i ((P(-, i), e.,)o e.,) . Then, using (fSSj) and taking 
into account that the operator (J-L — 7' id ) — ' acts continuously in L/2(F) and that the series in (|67p 
and (|68p converge in L^{F), we obtain the representations: 

S;(f)e, = («-7id)-' e-^^)(P(-,t), ej ) e,)dr. 

On the other hand, by the Elliptic Regularity Theorem and the Sobolev Embedding Theorem (see 
Section f3.1|) . the operator (H — jid)~ l acts continuously from L,2(F) into C k (F). Then we have 



= ||(^-7id)- i || B{L2 ,c7 fc )^. =i e- 2 ^(^o,e J )g < \\(H - 7id)- i || B(i2 ,^) ^Kllo, (70) 

llE^^^iHc* < ll(^-7idr Z || B(L2) ^)-|l / Yr =1 e ~ Hit ~ T) ( P ^ e ^e 3 dr\\ 
3 J $ 

< ||(^-7id)-'|| B(i2 ^) rilE^^^^^O^jOoe.Hodr 

< ||(ft - 7id)-'|| B(L2 ,c7*) jf ^E*!^^^-.*).^-)?^ 

< ||(^-7id)-'|| B(L2 ^) /V^-^llPMIIodT. (71) 

JO 

On the other hand, by Lemma 0n), 

||P(-, t) ||o < VVol(F, 9 ) ||(ft - 7 id )' (Q(-, t)) || c o < 0(1 + *) 2 'e 4A °* 
for some Q > 0. Then 

f e-^ {t - T) \\P(-,t)\\odT < Q f 6-^-^(1 + T ) 2l e iXoT dT 
Jo Jo 

<Qe-^(l + t) 21 / e^+ 4X ^ T dr = Q(l + t) 2l x\ " 1+4A ° ° T Ml ' (72) 

1 e-^H if 4A = -ah. 

From (f6T|) - ([72]) we get that claim (i) is valid. 

(ii) Recall that u = (u, e )o = «o + f °° qo(r) dr. From gU, jBTJ), ([gg ]> and (|7D|> - ([75 |1 we obtain 

n = e - A o*(u eo(x) + e(x,t)), 
where ||^(')i)||c7 fe = 0(e~ at ) for any € N, and a is the same as in claim (i). Hence 

Vu = e~ Xot (u Ve (x) + V0(s, t)) • 
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Furthermore, in view of and (|65|) . condition ([36]) yields uq 7^ 0. To show this, we calculate, 

using definition ([37]) of and condition ([36]) . 



o v&+ ~ \|/+ 

Using *ii- 3 e = (^/e^)(e /ii) 3 (eo/iio)'Uoeo and */e^ < e /u < 1/w-, e Q /u < l/u^, we get 

/ e Xot ( / ^(x)u- 3 (x,t)e (x)dx) dt < — (u , e )o / e A °*(u;_)" 3 (t) dt < (u , e ) . (73) 
Jo Jf Uq Jo 

Using w(x, t, t) := tu(x, t) + (1 — r)noe~ Ao< eo(x) = e _Aoi (uoeo(x) + r t)), we have 

— (V(\ogw(x,t,T)f) dr 

Z" 1 / V6(x,t) 6(x,t)(u Ve (x)+TV6(x,t)\ 
Jo {uoea(x)+r0(x,t) (u e (x) + t 9(x,t)) 2 J 

By the above, and the fact that 

inf {\u e (x) + T0(x,t)\ : x 6 F, t 6 [to, 00), r 6 [0,1]} > 
for to > large enough, follows claim (ii). □ 
Remark 2. For the extremal case of equality in (|36|) we have 

uo = Uq e {x), u = w (t)e (x), ^ = ^ + el(x). 
From ()3ip we obtain eo = const, hence /3 = const on F. Finally, u is a function of t. 



3.4 Nonlinear heat equation with parameter 

In this section the metric g, the connection V and the laplacian A smoothly depend on q which 
belongs to an open subset Q of R n . Consider the Cauchy's problem on a closed Riemannian 
manifold (F,g): 

d t u = Au + f(x,t,u,q), u(x,0,q) = u Q (x,q). (74) 

Here the function / is defined in the domain D = F x [0,T] x I x Q, where / C R is an interval, 
and Uo is defined in the domain D = F x Q and satisfies the condition: uq(x, q) € / for any x € F 
and q € Q. 

The next proposition is standard, for the convenience of a reader we prove it. 

Proposition 6. Suppose that f € C°°{D), uq € C°°{D), all partial derivatives of f and uq by 
x, u and q are bounded in D and D, and for any 5 £ Q there exists an unique solution u : 
F x [0, T] x Q — > R of the Cauchy's problem [7$ such that all its partial derivatives by x are 
bounded in F x [0, T] x Q. Then u(-,t, •) G C°°(F x Q) for any t € [0,T]. 

Proof. As is known, u(-,t,q) 6 C°°(F) for any q G Q, t G [0, T}. We should prove the smooth 
dependence on q of the solution u(x, t, q) and of all its partial derivatives by x for any fixed 
t € [0, T]. We shall divide the proof into several steps. 

Step 1: the continuous dependence of u(x,t,q) in q. To show this, take qo G Q and denote 
v(x, t, q) = u(x, t, q) — u(x, t, qo) and vo(x, q) = Uq(x, q) — uq(x, qo). Let us represent: 

f(x,t,u(x,t,q),q) - f(x,t,u(x,t,q ),q ) = F(x,t,q)v(x,t,q) + G(x,t,q) ■ (q - q ), 
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where 



fi 

F(x,t,q) = I d u f(x,t, u(x,t,q ) +rv(x,t,q), q + T(q- q ))dT, 



o 

G(x,t,q) = / V q f(x,t, u(x,t,q ) + Tv(x,t,q), q + r(q - q )) dr, 
Jo 

G (x,q) = / V q v (x, q + T(q- q ))dT, v (x,q) = G (x,q) ■ (q - q ). 
Jo 

Then the function v(x,t,q) is a solution of the Cauchy's problem: 

d t v = Av + F(x,t,q)v + G{x,t,q) ■ (q - q ), v\ t=0 = G {x,q) ■ (q - q ). (75) 
Then by the maximum principle of Proposition 

\v(x, t,q)\< w(t, q) V (x, t,q) € F x [0, T] x Q, (76) 
where w(t, q) is the solution of the following Cauchy's problem for the ODE: 

^ = F\w\ + G\q-q \, w(0, q) = G \q - q \ (77) 

with 

F= sup \F(x,t,q)\, G= sup \G(x,t,q)\, G = sup \G (x,q)\. 

(x,t,q)£Fx[0,T]xQ (x,t,q)sFx[0,T]xQ (x,q)&FxQ 

Then we get from (|76|) and (|77p the estimate 

\v(x, t,q)\< (G e Pt + (e Pt - 1)(G/F)) \q - q \, (x, t,q) 6 F x [0, T] x Q, 

which implies the claim of Step 1. 

Step 2: all the partial derivatives of u(x, t, q) by x depend continuously on q. Differentiating 
subsequently by x both sides of the equation and of the initial condition in (|74|) . we have the 
following Cauchy's problems for p a = d"u (a is the multi-index): 

d t pi = Api + d u f(x,t,u(x,t,q),q)pi+dlf(x,t,u(x,t,q),q), p i \ t=0 = d l x u {x, q), (78) 
d t pij = Apij + d u f(x, t, u(x, t, q),q)pij + (%f(x, t, u(x, t, q), q)piPj 

+ d u d l x f(x,t,u(x,t,q),q)pj + d^ 3 f{x,t,u{x,t,q),q), p itj \ t=0 = d x J u (x, q) (79) 

and so on. Applying the claim of Step 1 to these Cauchy's problems, we prove the claim of Step 2. 
Step 3: u(x, t, q) is smooth with respect to q. Take qo € Q and consider the divided difference 

x I * \ u(x,t,q + sy)-u(x,t,qo) 

dy(x,t,s) = (y el , s £ R). 

Denote 5^(x, s) = u °( x ' go+sy J "0(^.90) ^ j n ^ e same manner as in the proof of Step 1, we obtain the 
following Cauchy's problem for S y (x,t,s): 

d t 5 y = A5y + F(x,t,s)5 y + G{x,t,s)-y, 5 y \ t=0 = 5° y (x, s), (80) 

where 

F(x,t,s) = / d u f(x,t,u(x,t,q Q ) + r(u(x,t,q + sy) - u(x,t,q Q ))), q + sy) dr, 
Jo 

G(x,t,s) = / V q f(x,t,u(x,t,q ) +r(u(x,t,q + sy) - u(x,t,q ))), q + sy)dr. 
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Applying to the Cauchy's problem ([50]) the claim of Step 1, we conclude that S y (x, t, s) is continuous 
by s at the point s = 0, that is there exists the directional derivative 

d y (x,t,q ) = V q u(x,t,q ) ■ y = lim s ^ 5 y (x,t, s). 

Moreover, we see that d y (x,t, q) is the solution of the Cauchy's problem 

d t d y = Ad y + d u f(x, t, u(x, t, q), q)d y + V g f(x, t, u(x, t, q) • y, d y \ t=0 = V q u (x, q) ■ y. (81) 

Applying again the claim of Step 1 to this Cauchy's problem, we find that d y (x, t, q) = V q u(x, t, q)-y 
continuously depends on q for any direction y € M. n . This means that u(x,t,q) is C 1 -regular in q. 
Applying the above arguments to the Cauchy's problem ([5T]) . we obtain that u(x,t,q) belongs to 
C 2 with respect to q. Finally, we prove by induction that u{x,t, q) is smooth in q. 

Step 4. Applying all the arguments of Step 3 to the Cauchy's problems (|78p . (|79|) and so on, 
we prove that all derivatives of u(x, t, q) in x smoothly depend on q. □ 
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